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FOREWORD 


The  research  reported  herein  was  conducted  by  the  utaff  of 
the  Monsanto/Washington  University  Association  i.nder  tht  sponsor¬ 
ship  of  the  Advanced  Research  Projects  Agency,  Department  of 
Defense,  through  a  contract  with  the  Office  of  Naval  Research, 
N0001 4-67-C-0218  (formerly  N00014-66-C-0045),  ARPA  Order 
No.  873,  ONR  contract  authority  NR  356-484/4-13-66,  entitled 
"Development  of  High  Performance  Composites.  " 

The  prime  contractor  is  Monsanto  Research  Corporation. 

The  Program  Manager  is  Dr.  Rolf  Buchdahl  (phone  314-694-4721). 
The  contract  is  funded  for  $5,000,000  and  expires  30  April 

1970. 


ABSTRACT 


An  exploratory  investigation  of  mathematical  models  for  the 
description  of  wave  propagation  through  unidirectional  fiber- reinforced 
composites  is  reported.  A  self-consistent  model  using  three  concentric 
cylinders  with  different  materials  is  studied.  A  pure  shear  wave  propaga¬ 
tion  is  presented.  As  another  scheme  of  solving  the  dynamic  behavior, 
fiber  and  matrix  properties  may  be  simulated  by  periodic  functions, 
from  which  the  elasticity  theory  of  heterogeneous  media  can  be  applied. 


Introduction 


Dispersion  in  <m  elastic,  non-dissipstive,  wave-supporting  medium  is  a 
phenomenon  resulting  entirely  from  the  interaction  of  stress  waves  with  the  bound 
aries  of  the  medium.  Single  material  waveguides  admit  wave  reflections  at  their 
free  boundaries,  whereas  stratified  or  composite  waveguides  give  rise  not  only  to 
reflections  at  their  free  boundaries  but  also  to  reflections  and  refractions  of  waves 
at  interface  boundaries.  Mathematical  representation  of  dispers’-n  is  accomplished 
through  the  so-called  frequency  equation  which  relates  the  phase  velocity  of  the 
wave  to  its  wavelength.  When  there  is  no  dispersion,  the  phase  velocity  is  constant 
and  does  not  change  with  wavelength.  This  occurs,  for  example,  in  unbounded  elas¬ 
tic  homogeneous  media  where  waves  travel  with  the  dilatation  velocity 

1  I 

and/or  the  shear  wave  velocity  (Cs  =  (p/p)2).  Unbounded 

elastic  composite  media,  however,  cannot  support  waves  without  dispersion  because 
of  their  multiple  interface  boundaries.  Hence,  an  anisotropic,  homogeneous  repre¬ 
sentation  of  an  elastic  composite  medium  is  insufficient  for  predicting  its  response  to 
the  propagation  of  stress  waves.  This  illustrates  a  fundamental  difference  between  the 
static  and  dynamic  problem  for  composites. 

Presented  here  are  two  possible  modes  for  attempting  to  account  for  the  dis¬ 
persive  nature  of  composire  waveguides.  A  "self-consistent"  model  of  a  circular 
cylindrical  composite  waveguide  with  transverse  isotropy  is  proposed.  Three 
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concentric  circular  co-axial  cylinder*  constitute  the  model.  All  three  layer*  are 
elastic.  Isotropic  and  homogeneous  The  properties  of  the  outer  shell  are  the 
equivalent  properties"  of  the  composite  waveguide  being  modeled.  The  core  of 
the  model  has  the  properties  of  the  encased  fibers  of  the  composite  rod  and  the 
Intermediate  layer  has  the  properties  of  the  matrix  material  of  the  composite  rod. 

The  core  and  intermediate  layer  are  included  in  the  model  so  as  to  maintain  the 
identity  of  the  constituents  of  the  prototype  and  their  dimensions  can  be  adjusted 
to  control  the  dispersion  of  the  composite  rod.  Another  possibility  for  dispersive 
control  is  to  exchange  the  properties  of  the  core  and  intermediate  layer  of  the  model. 

A  second  mode  of  dispersive  representation  is  presented,  based  on  the  pre¬ 
mise  that  material  properties  which  can  be  expressed  as  a  function  of  position  give 
rise  to  dispersion.  The  potential  of  such  a  scheme  is  illustrated  by  a  simple  example. 
Material  property  functions  for  a  square  and  hexagonal  array  of  fiber  encasements 
are  proposed. 

I.  Self-Consistent  Model 

The  concept  of  a  "self-consistent"  model  has  been  introduced  and  a  develop¬ 
ment  of  its  frequency  equation  follows.  Conditions  under  which  a  considerably  sim¬ 
plified  yet  exact  solution  exists  are  demonstrated.  Phase  velocities  which  are  equal 
to  the  conventional  shear  wave  velocities  of  the  constituent  materials  can  exist  under 
certain  conditions  imposed  by  the  boundary  conditions. 
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Three-layer  stratification  of  the  model  is  necessary  so  that  enough  geometric 
parameters  are  available  to  adjust  the  equivalent  properties  of  the  composite  rod 
and  its  dispersion,  independently. 

In  what  follows,  material  properties,  physical  properties  and  configuration 
parameters  are  subscripted  with  a  \  ,  indicating  reference  to  the  core  when  j  1  , 
to  the  intermediate  layer  when  \  =  2  ,  and  to  the  outer  layer  when  j  =  3. 


1.  Solution  for  the  Displacements 

The  equations  of  motion  of  three-dimensional,  linear  elasticity  referred  to 
cylindrical  coordinates  are 

d2u.  2p  dai 

pj  =  V  2>V  V  ’  T1  H9'  +  2Hi 


(1) 


\2 

o  w. 

i  at2 


=  (X  +  2p.) 


&A. 

" 


2*  d 

5T 


(r«g.) 


2p.  d«L> 


where  u. 

1 


v  and  w  are  displacements  in  the  r  ,  9  and  z  directions,  respoc- 
1  1 


tively,  and 
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The  torsion  problem  for  circular  cylindrical  rods  Is  characterized  by  the 
following  assumed  solution. 


u.  =  w.  =  0 
1  1 

(3) 

«nd  v.  (r,  z,  t)  =  V.(r)  exp  [I  (yz  +  pt)] 

where  y  =  2*/ wavelength  and  p  =  2ir  (frequency).  Note  that  v.  Is  independent 
of  0.  Solution  (3)  in  conjunction  with  Equations  (1)  and  (2)  leads  to  a  single  ordinary 
differential  equation  in  V^(r). 
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Equation  (4)  is  a  Bessel  equation  of  order  1  and  has  the  fol lowing  solution: 


V.  =  A.J.  (k.r)  +  B.Y.  (k  r)  (5) 

1  1^1  1  1  1 

Solutions  for  each  layer  of  the  stratified  model  follow  immediately. 

V]  =  A]  J,  (k,r)  (6) 

V2  =  A2J,  (k2r)  +  B^  (k2r)  (7) 

■+ . 

V3  =  A3J,  <k3r)  ♦  B3Y,  (k3r)  W 


In  order  for  the  displacements  along  the  axis  of  the  model  to  remain  finite,  B1  -  0. 


The  k.  's  become 
1 


In  this  model  the  core  material  represents  the  encased  fiber  material  of  the  prototype 

and  the  subscripts  f  and  m  refer  to  the  fiber  and  matrix  material,  respectively. 
Bars  over  the  properties  in  k3  refer  to  the  "equivalent  properties"  of  the  prototype. 
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2.  Frequency  Equation 


Through  the  application  of  certain  boundary  conditions  along  the  lateral 
boundaries  of  the  model,  the  frequency  equation  is  obtained.  These  boundary  condi¬ 
tions  are: 


at  r  =  a. 


V,  =  V„  and  a  Q1  =  a 


'1  2 


r01  r02 


at  r  =  b, 


V„  -  V,  and  c  =  a 


2  ’3 


r02  r03 


at  r  =  c, 


ar03  =  ° 


00) 

01) 

02) 


The  form  of  solution  (3)  renders  all  other  stress  components  zero.  Application  of  the 
boundary  conditions  yields  the  following  set  of  homogeneous  equations  in  the 


con¬ 


stants  A.  and  B.  . 
1  1 


Vl(kla)  -  A^flcja)  -  B^O^a)  =  0 

Mf 

—  A^aJ^a)  “  2J1(k1a)]  -  A^k^J^a)  -  2J] (k2a)] 

"B2Ck2aY0(k2c)  "  2Yi(k2a)]  =  0 

A2Jl  (k2b)  +  BjY^kjb)  -  A3J1(k3b)  -  BgY^kgb)  =0  (1 3) 

T  {A2Ck2bJ0(k2b)  '  2Ji(k2b)1  +  B2Ck2bY0(k2b)  -  2^^)]} 

“A3Ck3bJ0(k3b)  '  2Jl(k3b)]  ”  B3[k3bVk3b)  “  2Y,  (k^)]  =  0 
A3Ck3cJ0(k3c)  "  2Jl(k3c)]  +  B3Ck3cY0(k3c)  “  2Y1(k3c)3  =  0 
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In  order  for  a  non-trivial  solution  for  the  A.  's  and  's  to  exist,  the  determinant 
of  their  coefficients  must  vanish.  Expansion  of  this  singular  determinant  vields  the 
frequency  equation  for  torsional  waves  travelling  in  the  self-consistent  mode^.  For 
each  value  of  y  an  infinite  number  of  values  of  p  will  result,  the  first  value  of 
which  is  that  corresponding  to  the  first  mode  of  propagation.  Dispersion  curves, 
sensitive  to  changes  in  material  and  physical  parameters,  are  thus  generated. 

As  previously  mentioned,  the  self-consistent  model  contains  sufficient  para¬ 
meters  to  adjust  equivalent  properties  and  dispersion  independently.  For  the  torsion 
problem,  the  equivalent  properties  will  depend  on  the  properties  of  the  constituent 
materials  of  the  prototype  and  its  fiber  volume  fraction,  Since  the  fiber  volume 
fraction  must  be  maintained  Throughout  the  self-consistent  model,  this  will  impose 
a  restriction  on  the  ratio,  b/a  ;  that  is,  the  ratio  of  the  outer  radius  of  the  intermedi¬ 
ate  layer  of  the  model  to  the  radius  of  the  core  will  depend  on  the  fiber  volume 
fraction,  £.  However,  with  the  properties  of  the  outer  layer  of  the  model  depending 
on  the  fiber  volume  fraction,  £,  and  not  on  the  interface  radii  a,  b  or  c,  the 
ratio  cA  is  free  to  be  chosen.  Variation  of  the  cA  ratio  allows  control  over  the 
dispersion  of  the  self-consistent  model  without  altering  its  equivalent  properties. 

3.  Torsional  Mode  with  No  Dispersion 

Under  the  conditions  that  both  constituent  materials  have  the  same  shear  wave 
velocity  (i.e.,  —  =  —  =  £)  a  considerably  simplified  solution  results.  Equations 

Pf  Pm  P 

(4)  reduce  to  one  equation  with  its  solution  describing  the  displacements  anywhere  in 
a  self-consistent  model. 


8 


ck .'jvfr).  o 


2  r  dr 


04) 


where  k  = 


r 

nr 


1/2 


.  ,  £  _  ^f  _  £ 

and  where 

f  m  K 


Equation  (14)  has  a  solution  of  the  form 


V(r)  =  AJj(kr) 


05) 


subject  only  to  the  condition  that  0^=0  at  the  free  boundary  of  the  model.  This 
condition  is 

kcJg(kc)  -  2Jj(kc)  =0  (16) 

The  first  solution  of  this  frequency  equation  is  k  =  0.  Therefore,  the  phase  velocity, 

1/2 

which  is  given  by  p/y  ,  is  equal  to  (p/p) '  .  Equa’son  (15)  does  not  represent  a 
solution  for  k  =  0  and  hence  Equation  (14)  must  be  re-iolved.  Therefore,  the  new 
solution  of  Equation  (14)  with  k  =  0  is 

V(r)  =  Ar  +  Br"1  (17) 

Since  stresses  and  displacements  must  remain  finite  at  the  origin  of  r 

V(r)  =  Ar  (18) 

It  is  interesting  to  note  that  Equation  (18)  represents  a  solution  for  the  first  mode  of 
propagation  for  the  torsion  problem  of  any  composite  waveguide  provided  that  the 
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free  boundary  is  circular  and  that  all  constituents  have  the  same  shear  wav.  velocity. 
This  is  true  for  any  number,  shape  and  eccentricity  of  the  encasements.  Not.  further 
that  this  mod.  exists  without  dispersion;  that  is,  the  ph.s.  velocity  is  constant. 


4.  Special  Solutions 

Suppose  now  that  the  phase  velocity  of  the  self-consistent  model  is  equal  to 
aithe,  the  sh.ar  wave  velocity  of  the  core  o,  the  shea,  wave  velocity  of  the  inter¬ 


mediate  layer,  with  ^  ^  ■  Since  the  following  development  is  conceptually 
pf  pm 


the  same  for 


either  case,  it  will  be  assumed  that  the  phase  velocity  is  equal  to  the 


1/2 


I  •»  rnr«  (C  =  (— )  )•  Under  this  circumstance.  Equations 

shear  wave  velocity  of  the  core  ;  / 


(4)  yield  the  following  solutions  for  the  self-consistent  model: 


(19) 
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Application  of  the  boundary  condition.  (10),  (1 1 )  and  (12)  yields  the  following  set 
of  homogeneous  equations  In  the  A.  's  and  B  's 

1  i 

aAl  -  J^k'ajAj  -  YjOc'ajBj  =  0 

Ck2aJ0(k2a)  '  2Jl(k2a^  A2  +  Ck2aY0(k2a)  “  2Y](k^a)]  =  0 
JjfkjbMj  +  Y1  (k^b)  B2  -  Jj  (kjb)  Ag  -  Yj  (k^b)  B3  =  0 

T\  Ck2bJ0(k2b)  '  2Jl(k2b)3  a2  +  Ck2bY0<k2b)  ‘  B2  j 

-  [Ic’bJoOc’b)  -  2J1(k^b)]  A3  -  [k^bY0(k'b)  -  2Y1(k’b)]  B3  =  0 
Ck3cJ0(k3c)  “  2Jl(k3c)1  a3  +  Ck3cY0(k3c)  -  2Y1(k^c)]  B3  =  0 


Again,  for  a  non-trivial  solution  of  the  Af  and  B. to  exist,  the  determinant 
of  their  coefficients  must  vanish.  This  singular  determinant  impose,  a  restriction  on 
the  values  of  y  (or  wavelength)  for  which  this  mode  can  exist. 

For  waves  travelling  at  the  velocity  of  shear  waves  in  the  intermediate  layer 

\Cp  =(^  ^  ^  8  S5mUar  freatmenf  wjH  W®ld  equation  which  impose,  a  restric¬ 

tion  on  the  value,  of  wavelength  for  which  this  mode  can  exist.  Therefore,  with  the 

exception  of  the  conditions  of  Art.  (3),  where  £  =  ^  ,  waves  travelling  with  the 

f  m 

conventional  shear  wave  velocity  of  either  constituent  material  will  be  allowed  only 


Hr  H 

discrete  values  of  wavelength.  That  Is,  unless  the  first  mode  of  torsional 

Kf  m 

wave  propagation  in  the  self-consistent  model  will  not  exist  without  dispersion. 

II.  Material  Properties  Which  Vary  with  Position  in  the  Waveguide 

Another  possible  method  of  predicting  the  response  of  an  elastic  composite 
waveguide  to  the  propagation  of  stress  waves  is  that  of  expressing  Its  material  prop¬ 
erties  as  functions  of  position.  This  method  has  the  advantage  of  allowing  the  equa¬ 
tions  of  motion  for  a  homogeneous  waveguide  to  be  used,  while  the  property  variance 
accounts  for  the  dispersive  nature  in  the  constitution  of  the  waveguide.  For  example, 
the  theory  for  homogeneous  rods  predicts  that  the  first  torsional  mode  is  propagated 
with  a  phase  velocity  equal  to  (£)l/2  and  that  this  mode  exists  without  dispersion. 

The  theory  for  the  self-consistent  mcdel  predicts  that,  in  general,  the  first  torsional 
mode  exists  with  dispersion.  A  similar  result  is  expected  for  the  theory  for  homo¬ 
geneous  waveguides  with  material  properties  which  vary  with  position. 

1.  One- Dimensional  Wave  Equation 

In  order  to  demonstrate  the  potential  of  the  material  property  variance  scheme 
consider  the  following  simplified  example.  Let  it  be  required  to  find  the  response 
of  plane  compressional  waves.  Let  the  waveguide  have  properties  which  somehow 
change  with  distance  in  the  direction  of  its  axis. 
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For  simplicity,  suppose  that  p  (mass  density)  is  constant  and  that  E  (modulus  of 
elasticity)  is  equal  to  E(z).  Then,  the  equation  of  motion  for  an  element  in  the 
waveguide  is 


p8z A 


d2u(z,  t) 


5zA 


daz(z,  t) 

5z 


0) 


where  A  is  the  cross-sectional  area. 


For  the  one-dimensional  problem,  Hooke's  Law  is 


•,w  = 


Equation  (1 )  then  becomes 


,  0  =  dE(z)  du(z,  t)  +  a2u(i,  t) 

4t2  *  41  Jz2 


Let  a  prime  denote  differentiation  with  respect  to  z  .  Then 


p - y{zt  0  =  E'(z)u'(z,  t)  +  E(z)u''(z,  t)  (2) 

dt 


Assume  a  solution  in  the  form 


u(z,  t)  =  U(z)  cos(wt) 


(3) 
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Then  Equation  (2)  reduces  to 

2 

d2U(z)  .  E£)  dUfei  +  «£  U(z)  -  0  <4> 

— r  +  ira  *  ^ 

dz 

.  ,h.  |atara|  dimensions  o(  the  waveguide. 

When  the  wavelength  is  large  compared  to 

approximate  solution  can  be  obtained  by  assuming  that 

.  \  (**) 
U(z)  =  A  cos(yz) 

where  y  =  2rr/wavelength. 

„„„  »,  V  .  —  -  W—  »>-"  “  -  *“  “ 

?  E'(z)A  /2itzv  +  ^ 

2  =  tjzjn  tan  (  )  + 

p  2irp  A 

Wore,  materia,  property . .  does  U-  to  dispersion  in  the 

Wev.  problem,  since, he  Phase  velocity  .Pen*  on  waueleng*.-.^-. 

E(I),  is  constant,  Ration  (6,  recces  to  the  conventional,  one-dimens.ona,  ph, 
velocity  for  compressional  waves,  Cp  fc/P  ) 

,  ,  ,  „,hi,  would  depend  on  the  form  of  E(z)  ),  «  does 

velocities  for  certain  values  of  z  th.s  w  . 

. . material  property  variance  can  l.ab  to  dispe.s, on  ,n  pro 

lems  which  *  no,  predict  dispersion  when  the  properties  constant. 

2^Variar«e_Function_for_S2uarep_Arra^_Fiber_Packin^ 

_ _ _  -  —  »•  —  ™  '  ““ 
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itself  to  a  relatively  simple  property  variance  function  and  will  be  treated  first.  The 
fiber-reinforced  waveguide  which  inspires  this  analysis  is  taken  to  be  a  cylindrical 
rod  with  its  fibers  parallel  to  the  longitudinal  axis  of  the  rod.  In  the  interest  of 
simplicity,  the  following  development  is  based  on  the  x,  y,  z  orthogonal  coordin¬ 
ate  system  with  z  co-axial  with  the  axis  of  the  waveguide.  If  A  (x,  y)  is  any  prop¬ 
erty  of  the  waveguide  which  varies  with  position  in  the  cross-section,  the  following 
property  variance  function  is  assumed  for  the  square  array. 

A(x,  y)  =  4  Ccos(— )  +  cos(^)]  +  B  (7) 

4  A  0 

X(x,  y)  is  a  continuous  function  so  that  for  any  point  in  the  cross-section  both  con¬ 
stituent  materials  contribute  to  the  property  in  a  meaningful  way.  Fiber  spacing  is 
the  same  in  the  x  and  y  directions  (for  simplicity)  and  the  property  function  has 
a  period  equal  to  this  spacing  in  both  directions.  The  constants  A  and  B  can  be 
determined  by  enforcing  certain  predetermined  behavior  of  A(x,  y).  These  condi¬ 
tions  are 

a.  X(ma,  na)  =  A^  (property  of  fiber)  at  the  center  of  each  fiber.  Therefore 

A  +  B  =  Af  (8) 

b.  A(ma, -y)  =  A(^?p,na)  =  A  at  points  intermediate  between  two 
adjacent  fibers,  where  A  is  the  conventional  equivalent  property  (for 
example,  the  law  of  mixtures).  Therefore 


B  =  A 


(9) 
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Equation  (9)  can  also  be  obtained  in  the  following  manner.  Consider  the  equation 

^//*X(x,y)dA  =  X  00) 

A 

This  equation  states  that  the  integral  or  X(x,  y)  over  the  cross-section  of  the  wave¬ 
guide,  divided  by  its  cross-sectional  area  is  equal  to  the  conventional  equivalent 
property  of  the  composite  configuration.  After  inserting  Equation  (7)  into  Equation 
(10)  and  integrating  over  the  limits  x,  y=  -p  to  +  p,  the  following  result  is  obtained 


If  the  dimensions  of  the  waveguide  are  an  integral  number  of  fiber  spacings,  then 
Equation  (11)  reduces  to 

X  =  B 

If  not,  however  the  fiber  spacing  is  small  compared  to  the  lateral  dimensions  of  the 
waveguide,  then  X  =  B.  Figure  2  shows  the  square  array  and  the  performance  of 
X(x,  y)  for  one  cycle  along  x  and  y  (a),  and  also  along  x  =  y  (b),  for  the  hypotheti¬ 
cal  case  when  Xf  =  10  psi ,  Xm  =  5psi,  £  =  0.8  and  X=Xf£  +  Xjl  -  f). 

When  the  law  of  mixtures  is  used  as  the  equivalent  property,  there  is  a  reason 
able  restriction  imposed  on  the  fiber  fraction,  £ .  At  points  equidistant  from  four  fiber 
inclusions,  X(ma/2,  na/2)  >  Xm  .  This  restricts  the  fiber  fraction  to  £  >  1/2.  If 
£=1/2,  X(ma/2,  na/2)  =  Xm  . 

3.  Variance  Function  for  Hexagonal  Array  Fiber  Packing 


A  more  complicated  property  variance 


function  for  hexagonal  packing  results, 
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as  might  be  expected.  Again,  using  X  as  a  general  property  which  varies  in  the 
cross-section  of  the  waveguide,  the  proposed  form  for  X(r,  0)  is 


X(r,  0)  =  A  cos  [3  +VT)  +  (3  -VT)  cos  60] } 

( 

+  |  C3  -  cos  603  [cos  1 E(3  +VT)  +  (3  -v^T)  cos  60] |  -  lj  +  C 


The  complication  of  this  function  requires  some  explanation.  In  order  that  X (r,  0) 
reduce  to  the  property  of  the  encased  fibers  at  the  center  of  each,  the  function  must 
have  a  period  which  varies  with  0  ;  hence  the  form  cos  (cos)  .  For  0=0,  v/3  , 
2k/3  ,  .  .  .  ,  the  period  is  equal  to  the  fiber  spacing,  a.  At  a/2  (intermediate  be¬ 
tween  two  adjacent  fibers),  X(r,  0)  =  X,  where  X  is  the  conventional  equivalent 
property  for  the  composite.  In  what  follows,  X  will  be  taken  to  be  the  law  of  mix¬ 
tures.  For  0=V 6,  3i*/6,  5t/6  ,  .  .  .,  X(r,  0)  has  a  period  aqua  I  to  VTa. 

Along  this  radius,  X(r,  0)  must  dip  to  a  minimum  at  r  =  ^^~a  (points  equidistant 
from  three  surrounding  fibers),  reach  a  secondary  maximum  at  r  =  ~^a  (points 

equidistant  between  two  adjacent  fibers),  dip  to  another  minimum  at  r  =  — V^~a 

3 

and  complete  its  cycle  with  a  maximum  of  X(r,  0)  =  Xf  at  r  =  VT a  .  This  behavior 
is  illustrated  in  Figure  4  for  the  same  hypothetical  case  used  for  the  square  array.  The 
sum  of  two  cos  functions  must  be  used  so  that  the  secondary  maximum  of  X (r,  0) 
can  be  obtained.  This  represents  one  of  the  conditions  for  the  determination  of  the 
constants  A,  B  and  C.  This  condition  can  be  separated  into  two  parts. 


a.  In  order  fcr  the  secondary  maximum  to  appear 


A  <  4B 


03) 
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>/ya 

b.  The  minimum  values  of  X(r,  0)  will  occur  at  r=  ^  and 

2  VT"  a  , 

r  =  — = when 


A  =  2B 


04) 


When  A  =  4B  ,  no  secondary  maximum  appears,  however  the  minimum  is  relatively 
flat  to  accomodate  fiber  spacing  large  compared  to  the  fiber  radiuf>,  b  .  Conditions 
(13)  and  (14)  lead  to  the  coefficient,  [3  -  cos  69]  . 


The  two  remaining  renditions  foi  the  determination  of  A,  B  and  C  are 

A  +  C  =  Xf  05) 

and 

-A  +  C  =  X  06) 

When  X  is  given  by  the  law  of  mixtures,  (X  =  X^£  +  Xm(l  -  £)  ),  conditions  (14), 
(15)  and  (16)  lead  to  the  following  expressions  for  A,  B  and  C. 

A  =  |[Xf(l  +  Xm(£-1)] 


B  =  ^[Xf(l  -  0  +  Xm(£-1)] 


07) 


C  =  |[Xf(5+4f)  -  4Xm(£-l)] 


HI.  ConcL.ions 

1)  Dispersion  is  predicted  for  the  first  torsional  mode  by  the  general  theory  for  the 
self-consistent  model. 
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2)  When  the  shear  wave  velocity  of  both  constituents  materials  is  the  same 


the  first  torsional  mode  for  the  self-consistent  model  is  propagated  without  dispersion 


3)  Particular  solutions,  with  discrete  allowable  values  of  wavelength  can  be  obtained 
when  the  phase  velocity  is  equal  to  the  shear  wave  velocity  of  either  constituent 


4)  The  self-consistent  model  appears  plausible  for  predicting  the  response  of  composite 
waveguides  with  transverse  isotropy  in  view  of  the  fact  that  its  equivalent  properties 
and  amount  of  dispersion  can  be  adjusted  independently. 


5)  Material  property  variance  functions,  when  used  in  wave  equations,  represent  a 
strong  possibility  for  predicting  the  response  of  composite  waveguides  to  the  propaga¬ 
tion  of  stress  waves. 


IV.  Brief  Outline  of  Immediate  Future  Work 

1)  Program  frequency  equation  for  self-consistent  model  and  obtain  dispersion  curves 
for  constant  equivalent  properties  and  various  c/b  ratios. 

2)  Develop  a  relation  between  c/b  and  the  amount  of  dispersion  in  a  prototype, 
composite  waveguide. 

3)  Proceed  to  solve  the  wave  equation  with  the  square  array  property  variance  func¬ 
tion.  More  sophisticated  functions  can  be  applied  later. 
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